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Abstract. A new method for constructing exact inhomogeneous universes is presented, that al- 
lows variation in 3 dimensions. The resulting spacetime may be statistically uniform on average, or 
have random, non-repeating variation. The construction utilises the Darmois junction conditions to 
v/^ ' join many different component spacetime regions. In the initial simple example given, the compo- 

^5 ■ nent parts are spatially flat and uniform, but much more general combinations should be possible. 

Further inhomogeneity may be added via Swiss cheese vacuoles and inhomogeneous metrics. This 
model is used to explore the proposal, that observers are located in bound, non-expanding regions. 
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p^ ■ while the universe is actually in the process of becoming void dominated, and thus its average ex- 

pansion rate is increasing. The model confirms qualitatively that the faster expanding components 
come to dominate the average, and that inhomogeneity results in average parameters which evolve 
differently from those of any one component, but more realistic modelling of the effect will need 

'^Vh ■ this construction to be generalised. 
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1 Motivation 

Up to now, the options for constructing exact inhomogeneous cosmologies were: (i) a small range of 
non-homogeneous metrics such as the Lemaitre-Tolman (LT) metric [1, 2], the Szekeres(S) metric 
[3, 4]/ their generalisations to non-zero pressure, the Lemaitre and Szafron metrics [1, 10], and 
a number of others with a less believable equation of state (EoS); (ii) the swiss cheese (SC) con- 
struction that inserts spherical structures into a Priedmann-Lemaitre-Robertson- Walker (FLRW) 
'background', including multi-level swiss cheese structures; etc.^ We here present a new way to 
construct exact inhomogeneous cosmologies that are arbitrarily inhomogenous out to all distances. 
They don't have an all-enveloping 'background' metric, or even an asymptotic one, and they could 
be inhomogeneous on any scale, but they can also be made to have the same 'average' everywhere. 
The scale and strength of the inhomogeneities can vary across the spacetime, or be kept statistically 
similar in all regions, and repeating patterns are also possible. It should be feasible to generalise 
this construction to other metrics than those considered here. 

In several recent works, Wiltshire [18-23] has discussed a range of unresolved questions about 
the assumptions underlying the standard approach to cosmological model building. In attempting 
to address them, he has propounded some alternative approaches, including the 'Cosmological 
Equivalence Principle', and some deep questioning of what we mean by 'averaging', and whether 
there is a well-defined relationship between an 'average' model and real observations. A key concept, 
that we attempt to model here, is the idea that void regions expand faster than cluster regions, 
and that, as time goes by, they occupy an increasing fraction of space, so the 'average' expansion 
rate becomes more and more dominated by the void expansion rate, while observers inhabit regions 
with little or no expansion, thus generating an apparent acceleration [24-27]. The model below is 



■"^See [5] for a survey of inhomogeneous cosmologies, [6, 7] for an overview and review of recent developments, and 
[8] for a well illustrated quick introduction and discussion of selected recent results. Also see the CQG issue [9] on 
inhomogeneous cosmologies. 

^For a good list of references see Grenon & Lake [11], footnotes 1 to 4. 

^The interesting method of Lindquist & Wheeler [12], has surface layers between the Schwarzschild cells, the 
nature of which is not clear. A recent non-vaccum generalisation [13, 14] is an approximate treatment. In the best 
current exact multi-black hole metric, there are multiple Reissner-Nordstrom bodies, which have gravitational and 
electromagnetic 'forces' exactly balanced [15, 16], with a A acceleration added [17]. 



relatively simple, being a first attempt of its kind, and it does not capture as many aspects of this 
particular proposal as we hope will be possible with subsequent generalisations. 

1.1 Inhomogeneity in LemaTtre-Tolman, Szekeres and Swiss Cheese IVIodels 

The Lemaitre-Tolman (LT) metric is spherically symmetric, but radially inhomogeneous, describing 
a ball of dust particles in comoving coordinates, for which the density and the dynamics depend 
on both time and radius. It contains the dust FLRW and Schwarzschild-Kruskal- Szekeres metrics 
as special cases, and is well suited to describing a black hole in a cosmological background. It is 
an excellent first approximation for non-linear gravitational collapse. As a cosmological model, it 
is certainly good for putting exact inhomogeneities, with a variety of scales, in an asymptotically 
uniform spacetime, but strong spherical inhomogeneities at very large radii are not observed. If one 
thinks of an LT model as an angular average of an inhomogeneous cosmology that is homogeneous on 
a sufficiently large scale, then the inhomogeneity should die off with radial distance. Nevertheless, 
it has seen extensive use in studying the non- linear evolution of cosmic structures, and in offering 
an explanation of the dimming of the supernovae. An asymptotically inhomogeneous model has 
been proposed to solve the horizon problem. 

The Szekeres (Sz) metric is even more interesting, since it has no Killing vectors, and thus 
doesn't suffer from the drawbacks of a high degree of symmetry noted above in the LT models. 
It has been used to model voids next to clusters, and even triple structures, and the range of 
possible structures it can describe is still not known. Nearly all studies have looked at the quasi- 
spherical class, which may be thought of as an LT model in which each spherical shell has been 
displaced relative to the others. This produces a dipole effect in the shell separation and in the 
density distribution around each shell. Thus it may be that this 'radial-inhomogeneity-plus-varying- 
dipole' eventually becomes rather unrealistic at large enough radii. The little studied quasi-pseudo- 
spherical class can be thought of as an irregular stacking of hyperboloids (pseudo-spherical shells) 
that vary in density and evolution, each of which contains an underdensity, or overdensity, due to 
the 'pseudo-dipole'. The principle inhomogeneity runs from one side of the universe to the other, 
with an extra variation snaking through the middle. It does not allow arbitrary inhomogneity in 
arbirary directions. Therefore this latter class appears to offer the possibility of very interesting 
structures and cosmologies, but it hasn't really been explored. 

In Swiss cheese models, one starts with an FLRW spacetime, and then cuts spherical holes 
in it, which may be filled with some other spacetime metric. For the construction to be valid, 
the Darmois junction conditions must be satisfied on the boundary between the two spacetimes; 
each boundary being a timelike 3-surface — that is, the history of a spherical 2-surface. Thus 
the possible interiors are restricted to metrics such as Schwarzschild, Lemaitre-Tolman, Vaidya, 
Szekeres. Lemaitre and Szafron interiors, that have non-zero pressure, are also possible but have 
not to our knowledge been used. Very often, the matching also restricts the FLRW EoS to be 
that of dust. Thus there is an obvious 'average' FLRW model, and the behaviour of the average is 
known from the start. 

The new construction has features significantly different from each of these. Further, it may be 
combined with all of the above, thereby much expanding the range of possible exact inhomogeneous 
cosmological models. 



2 Assembling Multiple 'Voids' & 'Clusters' in a Single Exact Spacetime 




Figure 1. A sample from a m.ulti-com.ponent inhomogcncous universe at one moment in 
time, showing several 'repeats' of a 'block' of component cuboids. In this illustration, each 
block (4 X 4 X 4 of them shown) has 8 components; 1 of type I, 3 of type II, 3 of type III, 
and 1 of type IV, and their sizes in each block are actually different, so they are never 
exactly repeated. These components are marked on the block at the near corner — the 3rd 
type II component (Ily) is of course hidden. In the main example of the paper, each 
component is a distinct Kasner-type region. The different K components have different 
expansion behaviour, so the fraction of a block they occupy changes with time. One could 
construct blocks with 3x3x3 components, or 5 x 2 x 7, or even component sequences in 
each of the x, y Sz z directions that never repeat. What matters is that the two expansion 
rates in the surfaces where two components join are the same in both component regions. 
Note that even when the set of expansion rates in a block repeats, the block sizes do not 
have to. This is the case illustrated above. 

In this section we construct a truly inhomogeneous universe, that does not have an obvious 
"background" metric. We use the Darmois junction conditions in General Relativity (GR) to join 
many regions of different matter content and evolution type. If the Darmois conditions are obeyed, 
the result is an exact solution of the Einstein Field Equations (EFEs). 

The essential idea is to do a 3-d tesselation, to fill space with an array of spacetime regions 
that are properly matched together at their boundaries. It is the need to ensure a proper 3-surface 
matching at each surface where component regions join, that renders this otherwise simple idea 
distinctly non-trivial. We begin by considering a kind of irregular cubic lattice, in which the basic 
'unit' or 'building block' is composed of 8 'pieces', or 'component regions', each one a different 



Kasner-type (K) metric [28] — see Fig. 1. The hope is that this construction may be generahsable 
to a variety of other metrics. In fact, an inhomogenous cosmology consisting of multiple slabs of 
FLRW and Kasner has been constructed in the "cheese slice" models of [29, 30], but the model 
presented here is much more realistic in that it is inhomogeneous in all 3 spatial dimensions, instead 
of just 1.^ 

Therefore, at any given moment, the 3-d patchwork is like a cubic lattice, but each of the 8 
components of a block has different expansion behaviour, so the relative sizes of each component 
will evolve greatly. The example illustrated in Fig. 1, needs 4 types of region, types I - IV. We 
could try to use type I regions to model 'voids', type II 'walls', type III 'filaments', and type IV 
'super clusters' ('clusters' for short). Each type I region adjoins 6 type II regions across a timelike 
3-surface; each type IV region adjoins 6 type III regions; each type II region adjoins 4 type III 
regions and 2 type I; and each type III region adjoins 4 type II regions and 2 type IV. To make sure 
the result is a regular solution of the EFEs, we will have to apply the Darmois junction conditions 
to the loci where 2 or more of these regions meet; see §4 for the details. 

In the simplest case, there is only one kind of type II & type III region, though oriented 
differently, and the basic block is cuboidal and repeats indefinitely. However, as subsequent sections 
will show, this much regularity is not required by the construction. In the more general case, the 3 
type II &: 3 type III regions may all be different metrics, and the expansion rates and sizes of the 
blocks and their components need not repeat, but rather have a random distribution of parameters. 
See section 5. 

3 The Kasner-type Metric 

The Kasner metric [28], and its non- vacuum generalisation, is a spatially flat, anisotropic Bianchi 
model of type I. It has a different 'expansion law' in each of 3 perpendicular directions. The metric 
is 

The Einstein field equations (EFEs), and the expansion are given by 

a/3 + /37 + 7" .3 2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 

The pressures are all proportional to the density, but different in the 3 perpendicular directions. 
The unit vector u"" = 51 is geodesic and comoving with the matter. 

The Minkowski Case To make the Riemann tensor zero requires e.g. 

a = = /3 = 7; or a = = /3,7 = l. (3.7) 
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*Soine properties were investigated in [31]. Although [32] adjoined different Bianchi I cosmologies on an initial 
time slice, and mentioned the case of FLRW next to Kasner, it is evident that a full Darmois-type matching was not 
achieved, since "delta-function discontinuities in the Riemann tensor" are found in section VI - i.e "surface layers" 
or "shock fronts" developed. A matching of planar dust (Ellis) metrics to vacuum was considered in [33]. 



The Vacuum Case The requirement p = = px = Py = Pz leads to (a, /3,7) being an equally 
spaced triplet round the Kasner circle: 

a=ii^. (3.8) 

^^l + 2cos|+W3)_ (3.9) 

It is not possible for more than two of a, /3 & 7 to be the same, so vacuum Kasner cannot be 
isotropic, unless it is Minkowski. The only vacuum case with two of them the same is {a, f3, 7} = 
{2/3,2/3,-1/3}. 

The Isotropic Case Putting 7 = /3 = a gives the spatially flat FLRW models, 

rr n, 3a2 2a-3Q2 / 2 \ ,^ _, 

where S is the FLRW scale factor. 

The Zero Pressure Case The requirement ^ = Px = Py = Pz '^s satisfied by the Minkowski case, 
the set of vacuum cases, and the a = 2/3 FLRW case. 

3.1 Physicality Conditions 

Since we will be assembling many different K regions to construct an inhomogeneous universe, it is 
useful to check their physical behaviour. The condition for non-negative density is easy to satisfy, 

p>0 -^ a/3 + /37 + 7Q >0 , (3.12) 

and that for non-negative pressure, 

Px > ^ /3 + 7 - (/3^ + /37 + 7^) > , 

^ \ ((1 - /3) - V(l - /3)(1 + 3/3) ) < 7 < ^ ((1 - /3) + V(l - /3)(1 + 3/3) ) , 
- 1/3 < /3 < 1 [ - 1/3 < 7 < 1] , (3.13) 

Py > -^ 7 + a- (7^ + 70 + 0^) > , (3.14) 

Pz>^ -^ a + /3-(a2 + a/3 + /32) >0 , (3.15) 

which is a fattened region around the Kasner circle, with a 3-lobed shape, a bit like a cardamom 
pod, is also not hard to satisfy. For a non-relativistic gas, we expect, 

Px<^ -^ 3(/3 + 7) - 3(/32 + /37 + 7^) < a/3 + /37 + 7a 
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7<-f3-4a-/3+ J9 + 12a - 6/3 - 20a2 _ 4a/3 + /32 
6 \ 

1 



& 7 > - f 3 - 4a - /3 - v/9 + 12a - 6/3 - 20a2 _ 4a/3 + /32 j ^ (315) 

Py < - -^ 3(7 + a) - 3(7^ + 7a + a^) < a/3 + ;S7 + 7a 



3 



7<-f3-4/3-a + \/9 + 12/3 - 6a - 20/32 - 4/3a + a^ 
6 \ 

1 



& 7 > - f 3 - 4/3 - a - v^g + 12/3 - 6a - 20/32 _ 4/3a + a2 j , (3.17) 

P^ < ^ ^ 3(a + ;S) - 3(a2 + a/3 + /32) < a/3 + /37 + 7a 



-^ 7<3-3a-/3 ^ . (3.18) 

Similarly, to ensure the sound speed less is than light speed, p < p, we obtain a similar set of 
conditions. These are not always so easy to obey, because, in the transition from dust FLRW, 
(a,/3,7) = (2/3,2/3,2/3) to Minkowski vacuum, (a,/3,7) = (0,0,0), the matter becomes increas- 
ingly stiff. Thus in the general case we have an exotic fluid. Nevertheless, as a first simple example 
of this type of construction, the above behaviour is not too bad. In any case, these component 
regions may themselves be thought of as averages over a more complicated matter distribution, so 
the above is only an effective bulk EoS (equation of state) . The fact that the pressures are different 
in different directions is only to be expected in regions with "pancake" or "cigar" expansion. More 
general metrics for the component spacetimes that are stitched together in such a patchwork, will 
no doubt allow us to improve this aspect. 

3.2 Cosmological Units 

We choose geometric units such that G = 1 = c, and for the remaining freedom, we specify that 
1 time unit = 10 gigayears. We call these cosmological geometric units, and the values of the 
cosmological time, length and distance units are: 

1 ctu = 10 Gy 

1 clu = 3.066 Gpc 

1 emu = 6.409 x 10^2 M© 

In these units, 100 km/s/Mpc is very close to 1/ctu. 

4 Junction Conditions 

We now implement the Darmois [34] junction conditions. If they are satisfied, the combined space- 
time metric is C^ and piecewise C^, and it may be shown [35]^ that, due to the Israel identities 
[36], the conservation laws SI yG^" = are satisfied even through the C^ junctions. 

For the 3-surface x = X = const, we define the surface S and the surface coordinates ^* to be 

x^ = (i,X,y,z) , (4.1) 

e = {t,y,z). (4.2) 

Here, greek indices range to 3, and latin indices range 1 to 3, but note that, in the context of 
junction conditions, index indicates the direction orthogonal to S, and is not necessarily time. 
We then calculate, in order, the surface basis vectors, their derivatives, the surface normal, the 
intrinsic metric and the extrinsic curvature: 

dx^ A o\ 

i^^y = IVT = ^ , (4.3) 

^^ Voooi/ 

n^ = (0,r,0,0) , (4.5) 

h = 9^^ ie^r (ejr = diag(-l,t2^t2^) , (4.6) 

Kij = -n. (^^ + T\, {ciY {e,A = . (4.7) 



di'di- 



1 



'Although this paper is about conservation faihng at a signature change surface, the non-signature-changing case 
was done first. 



The two manifolds we wish to join are labehed "+" and "— ", and are both Kasner-type (K) 
spacetimes: 

y+: x^ = (t+,x+,y+,z+) , 5/.. = diag(-l,t^+,t^+,tX+) , (4.8) 

V~ : x^ = (t„,x„,y„,z_) , 5/.^ = diag(-l,t"-,t!-,tI-) . (4.9) 

If the two surfaces to be identified, S4. & S_, are located at x = X+ &: x = X_, and the mapping 
between them is 

t+=t.=e , (4.10) 

2/+ = y_ = C^ (4.11) 

z+ = z.=e . (4.12) 

then the Darmois conditions, which use square brackets to denote the jump in a quantity across 
the junction, consequently require 

[h] =h\t-%]\i = ^ . (4.13) 

[Ki3\=K,j\t-Kij\^ = ^- (4.14) 



For our matching, these lead to 



/3+ = /3_ , 7+ = 7- • (4.15) 



This is consistent with the Israel requirement [36], which means that, for an observer moving with 
a timelike 3-surface, the pressures must match, but the density does not have to. Therefore, two 
K spaces may be joined across any pair of constant x surfaces if the two expansion indices in the 
plane of the surfaces are the same on either side of the junction. The expansion index perpendicular 
to S may differ in V+ &: y_. The same goes for matching pairs of constant y or z surfaces. Thin 
boundaries where the tangential pressures are different on either side are of course unrealistic; at 
the atomic level, one would expect streaming of particles and photons to blur the boundaries. 
However, junction condition methods are understood to be useful approximations to transitions 
that happen within a relatively thin region, for which the exact smooth equations are unsolvable. 

Now the Darmois conditions allow the regular junction of two manifolds at a pair of identified 
3-surfaces, the 3-surfaces being unbounded. In our case, each matching surface is bounded by 
other matching surfaces, so the size of the two identified surfaces must also be the same. Where 
2 cuboidal regions meet, in our array of different K regions, they must have 2 of the coordinate 
dimensions, Ax, Ay &: Az, in common — see below. 

Furthermore, where 4 cuboids touch at a 2-surface (a 1-d curve existing through time), the 
corner conditions [37] should be checked.^ However, we expect no difficulty, as the spatial coordi- 
nates in the K metric are Cartesian. Similarly, at the worldlines of the vertices where 8 cuboids 
meet, an appropriate modification of those corner conditions should be applied. 

5 Multi-Kasner Universes 

We are now ready to patch together Kasner-type (K) cuboids. Using the matching rules just 
obtained, a general, repeating 2x2x2 model would be as in Table 1: 



^At the level of bulk fluid parameters, the Darmois matching is complete. If however one introduces a kinetic 
theory description of the matter — far too complex for present purposes — then there would be further conditions 
needed to match all the modes. 

''See [31] for an application. 



Region 


Region type 


(a,/?, 7) 


{Ax,Ay,Az) 


I 


'void' 


(o,6,c) 


{A,B,C) 


IIx 


X 'wall' 


{d,b,c) 


{D,B,C) 


Ily 


y 'wall' 


{a,e,c) 


iA,E,C) 


IIz 


z 'wall' 


ia,b,f) 


iA,B,F) 


IIIx 


y-z 'filament' 


{a,e,f) 


iA,E,F) 


Illy 


z-x 'filament' 


id,b,f) 


{D,B,F) 


IIIz 


x-y 'filament' 


{d,e,c) 


{D,E,C) 


IV 


'cluster' 


id,ej) 


{D,E,F) 



Table 1. A repeating 2x2x2 multi-Kasner model. The region numbers correspond to 
those marked in Fig. 1. The triplet (a,/3,7) is the set of exponents in the metric (3.1), and 
the particular values a, 6, c, d, e, / may be chosen freely, subject to appropriate physicality 
conditions, as in section 3.1. Similarly the triplet (Ax, Ay, Az) gives the three coordinate 
dimensions of a region, as at the end of section 4, and A, B, C, D, E, F may be chosen 
freely. 

For a simple model, one may set 



/ 



d , C = B = A and 



E = D . 



(5.1) 



Thus, although this is not necessary, one of the 8 components could easily be an isotropic (i.e. 
FLRW) region, and we could even have two distinct FLRW regions per block. 

In fact, the building blocks do not have to repeat exactly. Referring again to Fig. 1, as we 
follow a line of component regions in the x direction, say, the particular pairs (/3, 7) and (S, C) must 
be the same for every region along that line, but need not be the same as adjacent lines. Similarly, 
following a line in the y direction, the particular pairs (a, 7) and {A, C) must be the same all down 
that row; etc. for following a z line. Thus, in each of the x, y k z directions respectively, we may 
have non-repeating sequences of expansion rates and widths: (ai, Ai), (q2, A2), (03, A3), • • • ; 
(^1, Bi), (/32, B2), (/33, B3), •••; (71, Ci), (72, C2), (73, C3), •••. Each (a^, Ai) stays constant 
in a constant {y, z) surface, and so on. Thereby we may create an arbitrarily inhomogeneous 
universe. If desired, these parameter values may be chosen from some statistical distribution. 

Furthermore, if there are FLRW regions, each one of them may be filled with a variety of swiss- 
cheese inhomogeneities, or even multi-level swiss-cheese inhomogeneites, and exact inhomogeneous 
models, such as the Lemaitre-Tolman and Szekeres metrics, can be used in the SC constructions. 



6 Averaging 

We now consider how to choose an FLRW model that best approximates this universe in some 
kind of 'average' sense, and we look at the evolution of the 'average' Hubble, deceleration and EoS 
parameters. 

The simplest approach is to do a volume-weighted average. When the 3-spaces are flat, this 
is quite simple. For simplicity, we consider a repeating 2x2x2 block, so that there is a uniform 
large-scale average. The volume occupied by each of the 8 component regions is just 



and its derivatives are 

V = AxAyAz{a + f3 + j) t"+/^+7-i ^ 

V = Ax Ay Az (a + /? + 7) (a + /? + 7 - 1) T+Z^+t-^ 
Thus the volume fraction of any given region i is 



fi 



T.UV, 



(6.1) 

(6.2) 
(6.3) 

(6.4) 



the volume-averaged expansion rate is 



e = ^i=^ ' = 3H , (6.5) 



and the volume-averaged deceleration parameter is 

_ „ 3(E;.iyi)(ELi'4) 

If we try to fit an FLRW model 



(6.6) 



«p = ^^^^ = 3H\1 + n,) , ^p = (2^^ y + k) ^ ^2(1 _ 2^ + f^^) ^ (6.7) 

to this inhomogeneous cosmology, then we would adjust the EoS, p = wp, so as to reproduce the 
actual expansion, i.e. 

and setting k = (since the K component regions are all flat), we have the effective EoS parameter 

W=i^. (6.9) 

For the simple example (5.1) above, these would be, 

(f3a j.2a+d j.a+2d f3d\ 

SJat^^^ (6.11) 



{f + f^) {a{a - 1)1" + d{d - l)^^) _ 



and for a single K component they are 



(6.12) 



e=51±^, , = -1 + ^^. „=i-_^. (6.13) 

t a + p + j a + P + 'y 

7 Model Details and Results 

We now present 3 variations of a specific case, a repeating 2x2x2 block. The details of models 
(A) to (C) are given in tables 2 to 4. In each, the type I & type IV regions are specified, and the 
remaining regions are fixed by the junction conditions. The repeating pattern allows an average 
to be calculated that is globally uniform, and thus the evolution of the component parts can be 
compared with the average. If the comoving size of a component region is A x S x C in the x, y 
& z directions, then the x size of the region today (t ~ 1 ctu) is ~ j4 clu; so A = 0.02 corresponds 
to ~ 60Mpc — the scale of voids, walls, filaments, etc. B k, C would have a similar order of 
magnitude. 



Region Region type Expansion rate triplet a-/3-j 



I 

IIx 

Ily 

IIz 

IIIx 

Illy 

Iliz 

IV 



'void' 
y-z 'wall' 
z-x 'wall' 
x-y 'wall' 
X 'filament' 
y 'filament' 
z 'filament' 
'cluster' 



(1/2,1/2,1/2) 
(2/3,1/2,1/2) 
(1/2,2/3,1/2) 
(1/2,1/2,2/3) 
(1/2,2/3,2/3) 
(2/3,1/2,2/3) 
(2/3,2/3,1/2) 
(2/3,2/3,2/3) 



P 



3 

11 
T2F 

11 
T2F 

11 
T2F 

10 
10 
10 

4 



{Px,Py,Pz 



( 1 1 1 ^ 

ViF' iF' W^ 

( 1 5 5 

ViF' 36l^' 36F 

/ 5 1 5 
V36l^' IP'' 36F 

/ 5 5 1 

V36l^' 36F' IF 

\^' 36F' 36F'' 

VseF'^' 36F'' 

V36F' 36F'^/ 

(0,0,0) 



e q 



3_ 

2t 



5 
3t 


4 
5 


5 
3t 


4 
5 


5 
3t 


4 
5 


11 


7 


6t 


11 


11 


7 


6f 


11 


11 


7 


6f 


11 


2 

t 


1 
2 



Table 2. Model (A). Type I regions (voids) are radiation FLRW, type IV regions 
(clusters) are dust FLRW. 



Region 


Region type 


Expansion rate triplet a-f3-j 


P 


{Px,Py,Pz) 


e 


Q 


I 


'void' 


(0,0,0) 





(0,0,0) 





- 


IIx 


y-z 'wall' 


(2/3,0,0) 





fo A A) 


2 

at 


7 
2 


Ily 


z-x 'wall' 


(0,2/3,0) 





(gt^^'O'Qt^) 


2 

3t 


7 
2 


IIz 


x-y 'wall' 


(0,0,2/3) 





(gF'gF'O) 


2 

3t 


7 
2 


IIIx 


X 'filament' 


(0,2/3,2/3) 


4 
9i^ 


(0,97^,9!^) 


4 
3t 


5 
4 


Illy 


y 'filament' 


(2/3,0,2/3) 


4 
9i^ 


(942,0,9^) 


4 
3t 


5 
4 


IIIz 


z 'filament' 


(2/3,2/3,0) 


4 
9i^ 


(gF'Qt^'O) 


4 
3t 


5 
4 


IV 


'cluster' 


(2/3,2/3,2/3) 


4 
3i^ 


(0,0,0) 


2 


1 
2 



Table 3. Model (B). Type 1 regions (voids) are (0,0,0) Minkowski vacuum, type IV 
regions (clusters) are dust FLRW. 



Region Region type Expansion rate triplet a-/3-7 p 



iPx,Py,Pz) 



e q 



I 

IIx 

Ily 

IIz 

IIIx 

Illy 

IIIz 
IV 



'void' 
y-z 'wall' 
z-x 'wall' 
x-y 'wall' 
X 'filament' 
y 'filament' 
z 'filament' 
'cluster' 



(1/10,1/10,1/10) 
(2/3,1/10,1/10) 
(1/10,2/3,1/10) 
(1/10,1/10,2/3) 

(1/10,2/3,2/3) 
(2/3,1/10,2/3) 
(2/3,2/3,1/10) 
(2/3,2/3,2/3) 



3 

looF 

43 

300*2 

43 
300*2 

43 
300^ 

26 
45F 

26 

45F 

26 
45F 

4 
3F 



' 17 
^ 100*2 ■ 

' 221 
^ 900*2 ■ 

' 221 
^ 900*2 ■ 

' 17 

^TooF' gooF' gooF' 
221 221 \ 



17 
100*2 : 

221 

900*2 : 

17 
100*2 ■ 

221 



17 - 
100*2 ' 

17 ' 
100*2 ' 

221 ' 
900*2 , 

221 < 



(0, 

( 



gooF' gooF' 



221 



221 



( 221 221 r,\ 

'^gooF' gooF' '' 
(0,0,0) 



3 

10* 

13 

15* 

13 

15* 

13 
15* 

43 
30* 

43 
30* 

43 
30* 

2 
* 



32 
13 

32 
13 
32 
13 

47 
43 

47 
43 

47 
43 

1 
2 



Table 4. Model (C). Type 1 regions (voids) are (1/10, 1/10, 1/10) FLRW near-vacuum, 
type IV regions (clusters) are dust FLRW. 

The behaviour of (B) & (C) should be quite similar, and we particularly focus on (C), even 
though the EoS p = (17/3)p is unrealistic, as it allows the evolution of both FLRW regions, as well 
as the average, to appear together on the graphs. The evolution of model (C) is plotted in Fig. 2. 
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nig,„(i) 

Figure 2. Evolution of the Dust-Near- Vacuum Model, with 0.02^ = ABC = DBF. Top 
left: evolution of the volume fraction of the 4 component types [type I (FLRW dust) - 
dashed, type II - dotted, type III - sohd, type IV (FLRW near- vacuum) - dot-dashed]; for 
types II & III there are 3 regions each with the fraction shown. Top right: evolution of the 
expansion [dashed - FLRW dust, dot-dashed - FLRW near- vacuum, solid - average]. 
Middle left: comparison of expansion rates [dashed - 0^^,/O/, dot-dashed - Qav/^iv]- 
Middle right: evolution of the deceleration [dashed - FLRW dust, dot-dashed - FLRW 
near- vacuum, solid - average]. Bottom: evolution of the effective equation of state 
parameter [dashed - FLRW dust, dot-dashed - FLRW near- vacuum, solid - average]. The 
units are cosmological geometric units, as given in §3.2, for which 1 ctu = 10 Gy, so the 
age of the universe is around 0.14 on the logj^g^ axis. The ordinate variables are 
dimcnsionless except for 0, which has units of fractional volume increase per ctu, so 
logj^g = 1 indicates a ten-fold increase in volume in 10 Gy, and a Hubble rate of 
65 km/s/Mpc corresponds to log^o © = 0.3. 

Figure 2 shows the time evolution of the volume fraction of each component region, the 
expansion rates of the average compared with regions I & IV, the average deceleration, and the 



11 



effective EoS parameter. This model clearly demonstrates how inhomogeneity can cause average 
parameters to evolve differently from what might be expected in a totally smoothed out universe. 
Although the vacuum or near- vacuum regions could model voids, the spatial flatness means 
they don't (or hardly) expand, whereas low density regions that develop in an initially expanding 
universe would retain their initial expansion — or rather the galaxies at their borders would. For 
greater realism, we would need to give the empty (type I) regions negative curvature (Milne model), 
and the high-density (type IV) regions positive curvature. The type II & III regions would need 
to be able to join them up across suitable surfaces and satisfy Darmois. Thus it is important to 
generalise this construction method by finding matchings for other types of component region. 

8 Discussion and Conclusions 

A new method for constructing exact inliomogeneous cosmological models is presented. The result- 
ing cosmologies may have arbitrarily strong inhomogeneities on arbitrary scales, and yet they may 
be made very random, or exactly homogeneous on average. They are not based on a 'background' 
or 'enveloping' metric that is effectively the large-scale average. The method can easily be com- 
bined with the swiss-cheese approach and exact inhomogeneous metrics, and thereby it provides for 
a much wider range of interesting comological models, with more degrees of inhomogeneity, than 
hitherto. Though the class of example models explored here is relatively simple, it does illustrate 
some significant features, and points the way to more complex possibilities. 

The above calculations have clearly demonstrated how the deceleration and equation-of-state 
parameters of the best-fit FLRW average model evolve with time, while those in each component 
of the construction are constant; and this even though the construction is exactly spatially uniform 
above a certain scale. This behaviour is due to the component regions having different expansion 
behaviour. Our finding therefore lends strong support to the contention that inhomogeneity means 
different regions evolve at different rates, and the evolution of the 'average' model is not the just 
the average of those rates, but strongly depends on which regions dominate the volume. Since 
we live in a region with little or no expansion, the average Riemann and Ricci fields can be very 
different from those felt locally. 

With the model presented here, the range of possibilities is not as large as one might like, 
since the spatial curvature of all the component regions is zero. More realistically, one would expect 
regions of both positive and negative spatial curvature, and one would expect some regions to 
expand and recollapse while others are ever-expanding. In this case, we expect the vacuum regions 
would expand fastest, and effective acceleration would emerge with time. However, it is a distinctly 
more interesting challenge to see how a collection of such regions could be patched together, using 
the technique presented here. While more general Bianchi models [38-40] seem to be the obvious 
extension, application of the Darmois conditions between positive and negative curvature regions 
is likely to result in much trickier constraints. Similarly, the limited range of constructions for this 
zero curvature case means that it is not possible to find completely reasonable equations of state for 
all the regions. We expect generalisation of this initial model will allow for some very interesting 
and realistic possibilities. Thus one should regard this current model as more illustrative than 
physical. 
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